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Abstract 

Adaptive wave model for financial option pricing is proposed, as a high-complexity 
alternative to the standard Black-Scholes model. The new option-pricing model, repre- 
senting a controlled Brownian motion, includes two wave-type approaches: nonlinear and 
quantum, both based on (adaptive form of) the Schrodinger equation. The nonlinear ap- 
proach comes in two flavors: (i) for the case of constant volatility, it is defined by a single 
adaptive nonlinear Schrodinger (NLS) equation, while for the case of stochastic volatil- 
ity, it is defined by an adaptive Manakov system of two coupled NLS equations. The 
linear quantum approach is defined in terms of de Broglie's plane waves and free-particle 
Schrodinger equation. In this approach, financial variables have quantum-mechanical 
interpretation and satisfy the Heisenberg-type uncertainty relations. Both models are 
capable of successful fitting of the Black-Scholes data, as well as defining Greeks. 

Keywords: Black-Scholes option pricing, adaptive nonlinear Schrodinger equation, 
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1 Introduction 



Recall that the celebrated Black-Scholes partial differential equation (PDE) describes the 
time-evolution of the market value of a stock option [BIS]- Formally, for a function u = u(t, s) 
defined on the domain < s < oo, < i < T and describing the market value of a stock 
option with the stock (asset) price s, the Black-Scholes PDE can be written (using the 
physicist notation: d z u = du/dz) as a diffusion-type equation: 

d t u = --(as) 2 d ss u - rsd s u + ru, (1) 

where a > is the standard deviation, or volatility of s, r is the short-term prevailing 
continuously-compounded risk-free interest rate, and T > is the time to maturity of the 
stock option. In this formulation it is assumed that the underlying (typically the stock) 
follows a geometric Brownian motion with 'drift' \x and volatility a, given by the stochastic 
differential equation (SDE) |3j 

ds(t) = ns(t)dt + as(t)dW(t), (2) 

where W is the standard Wiener process. The Black-Scholes PDE ([I]) is usually derived from 
SDEs describing the geometric Brownian motion ([2]), with the stock-price solution given by: 

s(t) = s (0)e^-^ t+tTW ^. 

In mathematical finance, derivation is usually performed using ltd lemma [I] (assuming that 
the underlying asset obeys the ltd SDE), while in physics it is performed using Stratonovich 
interpretation [6] (assuming that the underlying asset obeys the Stratonovich SDE [8]). 

The Black-Sholes PDE (pQ) can be applied to a number of one-dimensional models of 
interpretations of prices given to u, e.g., puts or calls, and to s, e.g., stocks or futures, 
dividends, etc. The most important examples are European call and put options, defined by: 

ucallO, *) = s-V(di) e~ TS - kAf(d 2 ) e~ rT , (3) 
u Put (s, t) = kAf(-d 2 ) e~ rT - sJV(-di) e~ TS , (4) 

^ (A) =K i+erf (^))- 

, Mi) +^(>- ln(f)+T(r-«-£ 

d2 



o-VT ' * o-VT 

where erf(A) is the (real-valued) error function, k denotes the strike price and 5 represents the 
dividend yield. In addition, for each of the call and put options, there are five Greeks (see, 
e.g. [9l HQ]), or sensitivities, which are partial derivatives of the option-price with respect 
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to stock price (Delta), interest rate (Rho), volatility (Vega), elapsed time since entering into 
the option (Theta) , and the second partial derivative of the option-price with respect to the 
stock price (Gamma). 

Using the standard Kolmogorov probability approach, instead of the market value of an 
option given by the Black-Scholes equation ([I]), we could consider the corresponding prob- 
ability density function (PDF) given by the backward Fokker-Planck equation (see [7]). 
Alternatively, we can obtain the same PDF (for the market value of a stock option) , using the 
quantum-probability formalism [1 1|, I12j . as a solution to a time-dependent linear or nonlinear 
Schrddinger equation for the evolution of the complex-valued wave ^—function for which the 
absolute square, is the PDF. The adaptive nonlinear Schrddinger (NLS) equation was 
recently used in |10| as an approach to option price modelling, as briefly reviewed in this 
section. The new model, philosophically founded on adaptive markets hypothesis |13|. 
and Elliott wave market theory |151 [16] , as well as my own recent work on quantum congition 
[17^118]. describes adaptively controlled Brownian market behavior. This nonlinear approach 
to option price modelling is reviewed in the next section. Its important limiting case with 
low interest-rate reduces to the linear Schrodinger equation. This linear approach to option 
price modelling is elaborated in the subsequent section. 
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Figure 2: Fitting the Black-Scholes put option with /3(u>)-adaptive PDF of the shock-wave 
NLS ip 2 (s,t) solution (fTU|) . Notice the kink near s = 100. 



2 Nonlinear adaptive wave model for general option pricing 

2.1 Adaptive NLS model 

The adaptive, wave-form, nonlinear and stochastic option-pricing model with stock price 
s, volatility a and interest rate r is formally defined as a complex-valued, focusing (1+1)- 
NLS equation, defining the time-dependent option-price wave function ip = ip(s,t), whose 
absolute square ^(s, i) | 2 represents the probability density function (PDF) for the option 
price in terms of the stock price and time. In natural quantum units, this NLS equation 
reads: 

id t i; = ~ad ss ^ - m 2 i>, (i = >/=!), (5) 

where /3 = /3(r,w) denotes the adaptive market-heat potential (see [19]), so the term V(ip) = 
— j3\il)\ 2 represents the ip— dependent potential field. In the simplest nonadaptive scenario j3 
is equal to the interest rate r, while in the adaptive case it depends on the set of adjustable 
synaptic weights {wj} as: 

(3(r,w)=r^w\ert(^y (6) 

Physically, the NLS equation © describes a nonlinear wave (e.g. in Bose-Einstein conden- 
sates) defined by the complex-valued wave function ip(s, t) of real space and time parameters. 
In the present context, the space-like variable s denotes the stock (asset) price. 



4 



The NLS equation ([5]) has been exactly solved using the power series expansion method 
[20\ 121] of Jacobi elliptic functions [22]. Consider the ip— function describing a single plane 
wave, with the wave number k and circular frequency uj: 

ij>(s,t) = 0(£) e 1 ^-"*), with £ = s-akt and 0(0 G R. (7) 

Its substitution into the NLS equation ([5]) gives the nonlinear oscillator ODE: 

f(0 + [w-^ 2 l^)+^(0 = 0. (8) 

We can seek a solution 0(£) for ([8]) as a linear function |21j 

0(0 = oo + aisn(£), 

where sn(s) = sn(s,m) are Jacobi elliptic sine functions with elliptic modulus m G [0, 1], such 
that sn(s,0) = sin(s) and sn(s, 1) = tanh(s). The solution of ([Sj) was calculated in [TO] to be 



±m^l sn(£), for me [0,1]; and 



<j){£) = ±w tanh(^), for m = 1. 
This gives the exact periodic solution of ([5]) as [TO] 



^(s,t) = ±m^jj^sn{s-akt)e i[ks -^ t(1+m2+k2)] , for mG [0,1); (9) 

ifo(s,t) = ± A l^-r tanh(s - akt) ^-^(2+^)] for m=1 ( 10 ) 

V 

where ([9]) defines the general solution, while (fTOj) defines the envelope shock-wav^ (or, 'dark 
soliton') solution of the NLS equation ([5]). 

Alternatively, if we seek a solution </>(£) as a linear function of Jacobi elliptic cosine 
functions, such that cn(s,0) = cos(s) and cn(s, 1) = sech(s)H 

HO = ao + aicn(£), 



A shock wave is a type of fast-propagating nonlinear disturbance that carries energy and can propagate 
through a medium (or, field). It is characterized by an abrupt, nearly discontinuous change in the character- 
istics of the medium. The energy of a shock wave dissipates relatively quickly with distance and its entropy 
increases. On the other hand, a soliton is a self-reinforcing nonlinear solitary wave packet that maintains its 
shape while it travels at constant speed. It is caused by a cancelation of nonlinear and dispersive effects in 
the medium (or, field). 

2 A closely related solution of an anharmonic oscillator ODE: 

<^."(s) + <^(s) + 3 (s) = 
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then we get [TO] 



± 




(}{w) 



a 



(}{w) 



a 



sec: 



cn 



)] 



for m = 1, 



for m G [0, 1); 



(12) 



(11) 



where (fTTj) defines the general solution, while ([12]) defines the envelope solitary-wave (or, 
'bright soliton') solution of the NLS equation ©. 

In all four solution expressions ©, (fTOj) . ([TT]) and (fT2|) . the adaptive potential /3(u>) is 
yet to be calculated using either unsupervised Hebbian learning, or supervised Levenberg- 
Marquardt algorithm (see, e.g. [23j [21]). In this way, the NLS equation ([5]) becomes the 
quantum neural network (see [18j). Any kind of numerical analysis can be easily performed 
using above closed-form solutions ipi(s,t) (i = 1, ...,4) as initial conditions. 

The adaptive NLS-PDFs of the shock-wave type (HOP has been used in [TO] to fit the 
Black-Scholes call and put options (see Figures Q] and [2]) . Specifically, the adaptive heat 
potential ([6]) was combined with the spatial part of (fTOl) 



while parameter estimates where obtained using 100 iterations of the Levenb erg-Mar quardt 
algorithm. 

As can be seen from Figure ([2]) there is a kink near s = 100. This kink, which is a natural 
characteristic of the spatial shock- wave (|13p . can be smoothed out (Figure [3]) by taking the 
sum of the spatial parts of the shock- wave solution (|10p and the soliton solution (|12p as: 



The adaptive NLS-based Greeks (Delta, Rho, Vega, Theta and Gamma) have been defined 
in [TO] , as partial derivatives of the shock- wave solution (|10p . 




(13) 




(14) 



is given by 
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Figure 3: Smoothing out the kink in the put option fit, by combining the shock-wave solution 
with the soliton solution, as defined by (|14|) . 



2.2 Adaptive Manakov system 

Next, for the purpose of including a controlled stochastic volatility^ into the adaptive-NLS 
model ([5]), the full bidirectional quantum neural computation model [18] for option-price 
forecasting has been formulated in [10] as a self-organized system of two coupled self-focusing 
NLS equations: one defining the option-price wave function ip = ip(s, t) and the other defining 
the volatility wave function a = cr(s,t): 

Volatility NLS : id t a = -^d ss a - (3(r,w) (\a\ 2 + |Vf ) a, (15) 

Option price NLS : id t ip = -^d ss ifj - 0(r, w) (\a\ 2 + |V>| 2 ) tp. (16) 

In this coupled model, the cr-NLS (I15p governs the (s, t)— evolution of stochastic volatil- 
ity, which plays the role of a nonlinear coefficient in (|16|) : the ^-NLS f)16|) defines the 
(s, t)— evolution of option price, which plays the role of a nonlinear coefficient in (I15p . The 
purpose of this coupling is to generate a leverage effect, i.e. stock volatility is (negatively) 
correlated to stock returns^ (see, e.g. [22]). This bidirectional associative memory effectively 
performs quantum neural computation [18] . by giving a spatio-temporal and quantum gener- 
alization of Kosko's BAM family of neural networks [28 \ 129] . In addition, the shock- wave and 
solitary-wave nature of the coupled NLS equations may describe brain-like effects frequently 

3 Controlled stochastic volatility here represents volatility evolving in a stochastic manner but within the 
controlled boundaries. 

4 The hypothesis that financial leverage can explain the leverage effect was first discussed by F. Black [26] . 
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occurring in financial markets: volatility/price propagation, reflection and collision of shock 
and solitary waves (see [30]). 

The coupled NLS-system (fT5]) - (fT6|) . without an embedded w— learning (i.e., for constant 
j3 = r - the interest rate), actually defines the well-known Manakov system^ proven by S. 
Manakov in 1973 [31] to be completely integrable, by the existence of infinite number of 
involutive integrals of motion. It admits 'bright' and 'dark' soliton solutions. The simplest 
solution of (fT5]) - (fT6|) . the so-called Manakov bright 2-soliton, has the form resembling that 
of the sech-solution (JT2J) (see [Ml US S3 S3 ESI HE ID]), and is formally defined by: 

Aoi(s,t) = 26csech(26(s + 4at))e- 2i(2a2 * +as - 2i,2i) , (17) 

where ip so i(s,t) = ( , s ), c = (ci,C2) T is a unit vector such that Icil 2 + I Co 1 2 = 1. Real- 

valued parameters a and b are some simple functions of (a, (3, k), which can be determined by 
the Levenb erg-Mar quardt algorithm. I have argued in [10] that in some short-time financial 
situations, the adaptation effect on /3 can be neglected, so our option-pricing model (fl~5j) - (fl~6j) 
can be reduced to the Manakov 2-soliton model (fT7|) , as depicted and explained in Figure [H 




Figure 4: Hypothetical market scenario including sample PDFs for volatility |<r| 2 and l^l 2 of 
the Manakov 2-soliton (]17p . On the left, we observe the (s, t)— evolution of stochastic volatil- 
ity: we have a collision of two volatility component-solitons, Si(s,t) and S2(s,t), which join 
together into the resulting soliton 5*2 (s, t), annihilating the Si(s, t) component in the process. 
On the right, we observe the (s, t)— evolution of option price: we have a collision of two op- 
tion component-solitons, Si(s,t) and S2(s,t), which pass through each other without much 
change, except at the collision point. Due to symmetry of the Manakov system, volatility 
and option price can exchange their roles. 



5 Manakov system has been used to describe the interaction between wave packets in dispersive conservative 
media, and also the interaction between orthogonally polarized components in nonlinear optical fibres (see, 
e.g. [321 133] and references therein). 
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3 Quantum wave model for low interest-rate option pricing 



In the case of a low interest-rate r <C 1, we have f3(r) <C 1, so V(tp) —> 0, and therefore 
equation © can be approximated by a quantum-like option wave packet. It is defined by 
a continuous superposition of de Broglie's plane waves, 'physically' associated with a free 
quantum particle of unit mass. This linear wave packet, given by the time-dependent complex- 
valued wave function ip = ip(s,t), is a solution of the linear Schrodinger equation with zero 
potential energy, Hamiltonian operator H and volatility a playing the role similar to the 
Planck constant. This equation can be written as: 

2 

iadttp = Hip, where H = ——d ss . (18) 

Thus, we consider the ^—function describing a single de Broglie's plane wave, with the 
wave number k, linear momentum p = ok, wavelength = 2ir/k, angular frequency ujk = 
o~k 2 /2, and oscillation period T& = 2ir/ujk = 4n/ak 2 . It is defined by (compare with |41|, 1421 

□21) 

Ms,t) = Ae< ks ~^ = Ae i{ks -^ = Acos(ks - —t)+Aism(ks - ^-t), (19) 

where A is the amplitude of the wave, the angle (ks — w&i) = (ks — 2 | L t) represents the phase 
of the wave ipk with the phase velocity: Vk = Wk/k = ak/2. 

The space-time wave function ip(s,t) that satisfies the linear Schrodinger equation (|18p 
can be decomposed (using Fourier's separation of variables) into the spatial part <f)(s) and 
the temporal part e~ lult as: 

i/)(a,t) = cj)(s)e- iujt = <l>(s)e-* Et . 

The spatial part, representing stationary (or, amplitude) wave function, (f>(s) = Ae lks , satisfies 
the linear harmonic oscillator, which can be formulated in several equivalent forms: 

0" + fcV = O, 0"+(^)% = O, </>"+(^)% = 0, 0" + ^ = O . (20) 

Planck's energy quantum of the option wave ipf. is given by: E k = OLO k = ^(ak) 2 . 

From the plane-wave expressions (|19|) we have: if)j.(s,t) = Ae*( ps ~ Ekt ^ — for the wave 
going to the 'right' and ipk(s,t) = Ae~v ( - ps+Ekt ^ — for the wave going to the 'left'. 

The general solution to (|18p is formulated as a linear combination of de Broglie's option 
waves (I19p . comprising the option wave-packet: 

n 

ijj{s,t) = ^Ci^ ki {s,t), (withneN). (21) 

8=0 
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Its absolute square, \ip(s,t)\ 2 , represents the probability density function at a time t. 

The group velocity of an option wave-packet is given by: v g = dwk/dk. It is related to 
the phase velocity Vk of a plane wave as: v g = v k — Xkdvk/dXk- Closely related is the center 
of the option wave-packet (the point of maximum amplitude), given by: s = tdojk/dk. 

The following quantum-motivated assertions can be stated: 

1. Volatility a has dimension of financial action, or energy x time. 

2. The total energy E of an option wave-packet is (in the case of similar plane waves) given 
by Planck's superposition of the energies Ek of n individual waves: E = naujk = ^(o~k) 2 , 
where L = no~ denotes the angular momentum of the option wave-packet, representing 
the shift between its growth and decay, and vice versa. 

3. The average energy (E) of an option wave-packet is given by Boltzmann's partition 
function: 

nE h 

(F] En= nE k e-^ _ E k 

En=O e bT ebT ~ 1 

where b is the Boltzmann-like kinetic constant and T is the market temperature. 

4. The energy form of the Schrodinger equation (|18|) reads: Eip = iadtip. 

5. The eigenvalue equation for the Hamiltonian operator H is the stationary Schrodinger 
equation: 

2 

H<f>(s) = E<t>{s), or E<t>{s) = -y^^(s), 



which is just another form of the harmonic oscillator (J20J). It has oscillatory solutions 
of the form: 

<p E {s) = c^^ 2 ^ 3 + c 2 e-- Vm ^ s , 
called energy eigen-states with energies Ek and denoted by: H(J)e(s) = Ek4>E(s). 

The Black-Scholes put and call options have been fitted with the quantum PDFs (see 
Figures [5] and [6]) given by the absolute square of (f2T|) with n = 7 and n = 3, respectively. 
Using supervised Levenb erg-Mar quardt algorithm and Mathematica 7, the following coeffi- 
cients were obtained for the Black-Scholes put option: 

a* = -0.0031891, t* = -0.0031891, h = 2.62771, k 2 = 2.62777, k 3 = 2.65402, 

k 4 = 2.61118, k 5 = 2.64104, k 6 = 2.54737, k 7 = 2.62778, c x = 1.26632, c 2 = 1.26517, 

c 3 = 2.74379, c 4 = 1.35495, c 5 = 1.59586, c 6 = 0.263832, c 7 = 1.26779, 
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Figure 5: Fitting the Black-Scholes put option with the quantum PDF given by the absolute 
square of (|21|) with n = 7. 



with (tbs = -94.0705a*, t BS = -31.3568**. 



Using the same algorithm, the following coefficients were obtained for the Black-Scholes call 
option: 



a* = -11.9245, t* = -11.9245, h = 0.851858, k 2 = 0.832409, 
k 3 = 0.872061, ci = 2.9004, c 2 = 2.72592, c 3 = 2.93291, 
with (T BS ~ 0.0251583ct*, t = -0.00838609** . 



Now, given some initial option wave function, ip(s,0) = ipo(s), a solution to the initial- 
value problem for the linear Schrodinger equation (|18p is, in terms of the pair of Fourier 
transforms (J 7 , J 7-1 ), given by (see |42j ) 



i>(s,t)=F- 1 [e-^^o)] =T 



e 2 



(22) 



For example (see [12]), suppose we have an initial option wave- function at time t = 
given by the complex-valued Gaussian function: 



g— as 2 /2gi<rfcs 



where a is the width of the Gaussian, while p is the average momentum of the wave. Its 
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Fourier transform, ifjo(k) = T[ip(s,0)], is given by 



'a 

The solution at time t of the initial value problem is given by 

1 f + °° -n <rfe 2 ^ a(k-p) 2 

iP(s, t) = -j== / e 1 ^-— *) e~^r- dk, 
which, after some algebra becomes 

ex P (- as2 -^f+; p2f ) 

ijj(s,t) = , (with p = ak). 

V 1 + iat 

As a simpler examplejf] if we have an initial option wave-function given by the real-valued 



8 An example of a more general Gaussian wave-packet solution of (|18|l is given by: 



( ~\a(s- s ) 2 - ±plt + ipo(s- s ) 



1 + iai V 1 + iat 

where So,Po are initial stock-price and average momentum, while a is the width of the Gaussian. At time t = 
the 'particle' is at rest around s = 0, its average momentum po = 0. The wave function spreads with time 
while its maximum decreases and stays put at the origin. At time —t the wave packet is the complex-conjugate 
of the wave-packet at time t. 
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Gaussian function, 



-72 



the solution of (|18p is given by the complex-valued tp— function, 



i;(s,t) 



exp(- 



2(l+ii) - 



From (I22h it follows that a stationary option wave-packet is given by: 
1 



000 



J- 



e* ks ij)(k)dk, where = -^^(s)]- 



As |0(s)| 2 is the stationary stock PDF, we can calculate the expectation values of the stock 
and the wave number of the whole option wave-packet, consisting of n measured plane waves, 
as: 



+00 



s\6(s)\ 2 ds 



and 



(k) 



k\^(k)\ 2 dk. 



(23) 



The recordings of n individual option plane waves (|19p will be scattered around the mean 
values (|23p . The width of the distribution of the recorded s— and A:— values are uncertainties 
As and Ak, respectively. They satisfy the Heisenberg-type uncertainty relation: 

71 

As Ak > -, 
~ 2' 

which imply the similar relation for the total option energy and time: 

n 

AE At > -. 

~ 2 

Finally, Greeks for both put and call options are defined as the following partial deriva- 
tives of the option ■0— function PDF: 



Delta = d s \ip(s,t)\ 2 = 

2i E?=i cjkj e^( is "5 iCT M)Abs [£™ Cj e ^( is -^'<) 



Abs' 



Vega = a (T |V(s,t)| 2 = 

-it YTj=i ^kj 2 e^O-s^Abs \j2j =1 Cj e^-^Ayt) 
Theta = d t \4>(s,t)\ 2 = 

-ia YTj=i Cjkj 2 e^'( is -5 iCTA; ^)Abs \Y%=i c j M*"-*™'*) 



Abs' \jyi= 1 c j M is -^ t ) 



13 



Gamma = d ss \ip(s, t)\ 2 = 
E"=i c,-fcj e^( is ~^ i,jfc ^)) 2 Abs' [E"=i Cj e^C" - ^^*) 



Abs' 

2_ 

Abs" 



^ n =1 c i e^( is ~5 iCTfc j*) 



where Abs denotes the absolute value, while Abs' and Abs" denote its first and second 
derivatives. 



4 Conclusion 

I have proposed an adaptive-wave alternative to the standard Black-Scholes option pricing 
model. The new model, philosophically founded on adaptive markets hypothesis |13t [T3] 
and Elliott wave market theory [151 I16j . describes adaptively controlled Brownian market 
behavior. Two approaches have been proposed: (i) a nonlinear one based on the adaptive 
NLS (solved by means of Jacobi elliptic functions) and the adaptive Manakov system (of 
two coupled NLS equations); (ii) a linear quantum-mechanical one based on the free-particle 
Schrodinger equation and de Broglie's plane waves. For the purpose of fitting the Black- 
Scholes data, the Levenberg-Marquardt algorithm was used. 

The presented adaptive and quantum wave models are spatio-temporal dynamical systems 
of much higher complexity [25] then the Black-Scholes model. This makes the new wave 
models harder to analyze, but at the same time, their immense variety is potentially much 
closer to the real financial market complexity, especially at the time of financial crisis. 
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